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Abstract 

We extend the recent computation of Higgs boson production in association with a jet through 
next-to-next-to-leading order in perturbative QCD by including decays of the Higgs boson to 
electroweak vector bosons. This allows us to compute fiducial cross sections and kinematic distri¬ 
butions including realistic selection criteria for the Higgs boson decay products. As an illustration, 
we present results for pp^H + j^yy + j closely following the ATLAS 8 TeV analysis and for 
pp ^ H + j ^ W~^W~ + j ^ e'^pTvv + j va a. CMS-like 13 TeV setup. 
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I. INTRODUCTION 


Studies of the Higgs boson discovered by the ATLAS and CMS collaborations [1, 2 ] will 
be at the focus of the experimental program during Run II of the LHC. The interpretation 
of future measurements of Higgs boson production and decay rates in terms of Higgs boson 
couplings to matter and gauge helds and Higgs boson quantum numbers will rely on the 
comparison of measured event rates and kinematic distributions with results of theoretical 
modelling of such processes in the Standard Model. It is hoped that such comparisons will 
help to elucidate the nature of the Higgs particle and explore the mechanism of electroweak 
symmetry breaking in detail [3, 4]. 

Recently, the ATLAS collaboration made an important step forward in presenting the 
results of analysis of the Higgs boson production and decay at the LHC. Indeed, in contrast to 
many other Run I LHC measurements, the ATLAS collaboration measured hducial volume 
cross sections and a variety of kinematic distributions in the process pp ^ H + jets [5]. The 
Higgs boson was identihed through its decay to photons, H 77 . 

Fiducial volume measurements allow for a direct comparison between data and theoretical 
predictions minimizing extrapolation uncertainties. Although many of the ATLAS hducial 
volume measurements are currently limited by statistical uncertainties, this will certainly 
change in the current run of the LHC. Therefore, the important issue for the near future 
is the availability of highly accurate theoretical predictions that can be used to describe 
complicated hducial volume measurements. 

We will now summarize the most advanced hxed order computations related to Higgs 
boson production and decay at the LHC. The inclusive production cross section of the 
Higgs boson has recently been computed through next-to-next-to-next-to leading order in 
perturbative QCD [ 6 ]. This computation refers to the total cross sections and can not be 
used for the direct comparison with hducial volume measurements without extrapolation. 
The computation oi H + j production at the LHC has recently been extended to next-to- 
next-to leading order (NNLO) in perturbative QCD, in a fully diherential manner [7, 8 ]^. 
Unfortunately, in Refs. [7, 8 ] decays of the Higgs boson were not considered; for this reason 
the comparison of the results of Refs. [7, 8 ] with the results of the hducial measurements 

^ For earlier partial results, see [9, 10]. 
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is also not possible. When the Higgs boson is produced in association with two and three 
jets, the NLO QCD computations provide the state-of-the-art results [11-13]. In those NLO 
QCD computations decays of the Higgs bosons are routinely taken into account. We note 
that the NNLO QCD computations oi pp ^ H + j [7, 8 ] combine the NNLO prediction for 
the exclusive H + j cross section with the NLO QCD prediction for the exclusive H -|- 2j 
cross section and the LO prediction for the H + 3j cross section, making them particularly 
suitable for studying Higgs boson production in association with different number of jets in 
a consistent way. 

It is relatively straightforward to extend the fully-differential pp ^ H + j computation 
reported in Ref. [7] to include decays of the Higgs boson into electroweak gauge bosons since 
the Higgs boson is a spin-zero particle and no spin correlations need to be considered. This 
is what we do in this paper for a variety of the Higgs boson decay modes. Once this is done, 
it becomes possible to calculate fiducial volume cross sections and kinematic distributions 
and directly compare with experimental measurements. The very fact that it is possible 
to do that through next-to-next-to-leading order in the expansion in the strong coupling 
constant, represents an impressive milestone in an application of perturbative QCD to the 
description of hard collisions at the LHC. 

Our paper is organized as follows. In Section H, we briefly summarize the theoretical and 
experimental setup. In Section HI, we present the results for fiducial volume cross sections 
and kinematic distributions at the 8 TeV LHC for the 77 —?■ 77 decay mode and at the 
13 TeV LHC for the 77 —>■ WW* —)■ vv decay mode. We also compare the results 

of the hducial volume computation for the 77 hnal state with the results of the ATLAS 
measurement. We conclude in Section IV. 

II. THE SETUP 

A. Theory 

We begin by summarizing the theoretical framework that we use in the computation. We 
work in an effective held theory obtained by integrating out the top quark. We employ the 
method of improved sector decomposition developed in Refs. [14-16]. This method is based 
on the factorization of scattering amplitudes in soft and collinear limits and on particular 
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way of splitting the phase-space into sectors, where soft and collinear singularities are easily 
identihed. For this calculation, we require a large number of matrix elements that are used to 
construct differential cross sections. In particular, we need the two-loop virtual corrections to 
the partonic channels gg —?■ Hg and qg —?■ Hq-, the one-loop virtual corrections to gg —)■ Hgg, 
gg —)■ Hqq^ qg —)■ Hqg, qq —)■ HQQ and the double real emission processes gg —)■ Hggg, 
gg —j- Hgqq, qg —Hqgg and qg — HqQQ, where the QQ pair can be of any flavor. The 
helicity amplitudes for all of these processes are available in the literature. The two-loop 
amplitudes were computed in Ref. [17]. The one-loop corrections to the four-parton processes 
are known [18]. For hve-parton tree-level amplitudes, we use compact results obtained using 
BCFW recursions [19]. 

It is non-trivial to combine processes with different particle multiplicities as required 
for any NNLO QCD computation. Our method for doing that is described in [9]; we do 
not repeat that discussion here. As already mentioned in the Introduction, the inclusion of 
Higgs boson decays is straightforward since the Higgs boson is a spin-zero particle. The only 
technical issue that arises is a signihcantly larger phase space that needs to be considered and 
the ensuing difficulties with the Monte-Carlo integration. However, the numerical challenges 
that appear in fiducial volume computations turn out to be not prohibitive. In particular, 
we find that for the —)• 77 decay mode we need roughly the same amount of statistics as 
for the stable Higgs case, while for H ^ 41 the amount of statistics should be increased by 
a factor between 2 and 4. 

B. R — 7 - 77 

We continue by listing the selection criteria employed by the ATLAS collaboration [5] . We 
are interested in the process pp —)■ H+j, where the Higgs boson decays to two photons. Final 
state jets are defined using the anti-fcj_ algorithm [20] with AR = 0.4 and p±j > 30 GeV. 
Jets are required to have rapidities pj in an interval —4.4 < pj < 4.4. The two photons 
from the Higgs decay must have the transverse momenta > max (25 GeV, 0.35 
and p ±,72 > max (25 GeV, 0.25 respectively, where is the invariant mass of the 

two photons. In our calculation, the Higgs boson decays are described in the narrow width 
approximation, so that we always have = mn = 125 GeV. Then, the above conditions 
imply > 43.75 GeV and ^±,72 > 31.25 GeV. ATLAS requires that the two photons 
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are in the central region of the detector \y^\ < 2.37, but no photons are allowed to be 
in the rapidity interval 1.37 < \y^\ < 1.56. However, when presenting the results of the 
measurements [5], the ATLAS collaboration corrects for the second condition and, for this 
reason, we do not account for it in our calculation.^ It is required that photons and jets are 
well-separated AR^j > 0.4. Finally, we take the branching ratio for the Higgs boson decay 
to two photons to be 77 ) = 2.35 x 10“^. 

C. H ^W+W- ^ e+y-vi? 

As the second example, we consider H + j production at the 13 TeV LHC. The Higgs 
boson decays to yTvv hnal state through a pair of W bosons. To identify selection criteria, 
we apply kinematic cuts similar to those employed by the CMS collaboration in their studies 
of the H —)■ H^+H^“ production at the 8 TeV LHC [21]. We dehne jets using the anti-fc_L 
algorithm with AR = 0.4. Jets are required to have transverse momentum > 30 GeV 
and be in the rapidity interval —4.7 < yj < 4.7. The harder of the two charged leptons must 
have transverse momentum p±^i > 20 GeV; the softer one must have p±^i > 10 GeV. The 
transverse missing energy in the event should exceed i7_L,miss > 20 GeV. Other cuts that 
we employ are i) the cut on the dilepton invariant mass mu > 12 GeV; ii) the cut on the 
transverse momentum of the dilepton pair pj_ ii > 30 GeV and iii) the cut on the transverse 
mass of the two W bosons m± = \^/2p \ gE \ mi^A l - cos > 30 GeV. 

III. THE RESULTS 

A. 77—7-77 

We consider production of the Higgs boson in association with a jet at the 8 TeV LHG and 
compute the hducial volume cross section and the kinematic distributions using the ATLAS 
selection criteria described in Section H. We begin with the hducial volume cross section. 
For events that contain the Higgs boson and at least one jet, the ATLAS collaboration 

^ We note that it is unclear to us why correcting for the missing rapidity region is a worthwhile thing to 
do since such theoretical correction defies the original goal of comparing theoretical fiducial volume cross 
sections with the results of experimental measurements. 
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Figure 1: Left pane: Fiducial cross section for pp —)• + j at 8 TeV LHC as a function 

of the cut on the jet transverse momentum. The inset shows ratios of differential cross sections 
at different orders in perturbation theory for the factorization and the renormalization scales set 
to the mass of the Higgs boson. Right pane: fiducial cross sections for exclusive jet bins. The 
selection criteria are described in the text. 


obtains [5] 


TeV) = 21.5 ± 5.3(stat.) ±1'^ (syst.) ± 0.6(lumi) fb. (1) 


This result is significantly higher than what the hxed order computation predicts. For the 
8 TeV LHC we obtain inclusive fiducial H ( 77 ) + j production cross sections^ 


_ c ^q+2.32 ri fid _ -7 QC+l -''’6 fu ^fid _ q ^[-+0.58 rL 

'^LO ~ <^NLO ~ ' •y8_^ 40 ID, UnNLO ~ ^'^^-0.82 


( 2 ) 


where the central value corresponds to the factorization and renormalization scales set to a 
common value /i = mn and the upper (lower) values to = mH/‘^ (r = 2m h), respectively. 
Theoretical results indicate reasonable convergence of the perturbative expansion. Indeed, 
for /i = mu/2 the NLO cross section is higher than the LO one by 26%, and the NNLO 
cross section is higher than the NLO one by only 3%. The situation is worse but similar for 
pi = mn, where the NLO (NNLO) corrections amount to 47 (18)%, respectively. 

It is interesting to point out that the quality of the perturbative expansion for pp ^ H+j 
appears to be somewhat better than for the inclusive Higgs boson production. As was 
pointed out earlier in the context of NLO computations, this feature may be related to 
a cancellation of Sudakov logarithmic corrections 0{a'^ p±j/m^) that are present in 

^ We remind the reader that we work in the m* —>■ oo limit. Moreover, we do not include subleading qq 
channels beyond NLO, c/. Refs. [7, 8]. 
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the computation of the H + j cross section, and other sources of large corrections that 
contribute to the inclusive rate. It was argued that this cancellation is accidental and for 
this reason can not be expected to continue in higher orders of perturbation theory. Our 
calculation questions this assertion. Indeed, it clearly demonstrates that, through the NNLO 
in perturbative QCD, there is no indication that hxed order perturbation theory for H + j 
production with the cut on the jet transverse momentum > 30 GeV breaks down. 

As a further illustration of this point, we show in the left pane of Fig. 1 the cross sections 
for pp ^ H + j process as a function of the jet transverse momentum cut p_L,cut- It follows 
from this plot that both the NLO and the NNLO QCD corrections are moderate for all 
values of p_L,cut- In particular, the convergence of perturbative series for the cross section 
with p_L,cut = 30 GeV does not appear to be signihcantly worse than the convergence of 
perturbative predictions for larger values of the transverse momentum cut. 

It is interesting to point out that in Ref. [5] a summary of the theoretical results for 
H + j hducial volume cross sections was presented, based on exact NLO QCD computations 
and various re-summations of potentially enhanced terms [22-24] . Compared to these three 
results, our result is somewhat higher than the result of Ref. [22] and it is lower than the 
two results based on the calculations of Refs. [23, 24].^ On the other hand, the uncertainty 
of our result is signihcantly smaller. 

Nevertheless, the difference between the measured pp ^ H + j cross section by the AT¬ 
LAS collaboration and our theoretical prediction is striking ~ the central value measured 
by ATLAS is higher than our results by a factor 2.1 — 2.5, depending on the choice of the 
renormalization and factorization scale. If one accounts for the uncertainty of the experi¬ 
mental result, this difference translates to, approximately, 2.4 standard deviations and is not 
statistically signihcant. On the other hand, it is interesting to remark that the mismatch 
in the H + j channel is stronger than the mismatch in the inclusive pp ^ H cross section, 
where the experimental cross section exceeds the theoretical one by a factor of 1.4 [5]. It will 
be interesting to watch how results of these measurements will evolve in the future especially 
since, thanks to the availability of the NNLO QCD computation, the precision of theoretical 
calculations is quite high. 

We can recast our results in Eq. (2) into predictions for acceptances, at different orders 
^ Note however that finite mass effects can increase the cross section by about 6% [25, 26]. 


7 




Figure 2: Higgs boson rapidity (left) and transverse momentum (right) distributions at the 8 TeV 
LHC. The insets show ratios of differential cross sections at different orders in perturbation theory 
for the factorization and the renormalization scales set to the mass of the Higgs boson. 

in perturbation theory. An acceptance is defined as the ratio of a hducial to total cross 
section A = agd/o' for H + j production. When ratios of cross sections are computed, 
many sources of theoretical uncertainties cancel out and it is in general not possible to 
properly estimate the uncertainty of the result by changing factorization and renormalization 
scales within a prescribed interval. For this reason, it is useful to know several orders in 
the perturbative expansion of the acceptance, to estimate the precision with which it can 
actually be predicted. For the 8 TeV LHC and the ATLAS setup, we hnd 

Alo = 0.594(4), Anlo = 0.614(3), Annlo = 0.614(4). (3) 

The perturbative expansion for the acceptances exhibits good convergence. Indeed, by com¬ 
paring the central values, we hnd that the NLO acceptance is larger than the LO acceptance 
by 3 percent, whereas there is no change going from NLO to NNLO. 

Another interesting quantity is the exclusive cross section for hxed number of jets. The 
corresponding results are shown in the right pane of Fig. 1. We observe good convergence 
of the perturbative expansion for the exclusive H + j and H -|- 2j production cross sections 
at 8 TeV LHC. We can not discuss the perturbative behavior of the H + 3j cross section 
since it enters our computation only at leading order in perturbative QCD. 

We now turn to kinematic distributions studied by the ATLAS collaboration. They can 
be divided into three categories: transverse momentum and rapidity distributions of the 
Higgs boson; transverse momentum, rapidity and the transverse energy distributions of the 
accompanying QCD radiation and, hnally, kinematic distributions of individual photons. 




























Figure 3: Rapidity and transverse momentum distributions of the most energetic jet at the 8 TeV 
LHC. The insets show ratios of differential cross sections at different orders in perturbation theory 
for the factorization and the renormalization scales set to the mass of the Higgs boson. 

The latter includes the transverse momentum and the rapidity distributions as well as the 
distribution of the photon decay angle in the Collins-Soper reference frame. We can compute 
all these kinematic distributions through NNLO in perturbative QCD, using exactly the same 
setup that the ATLAS collaboration employs in the actual measurement. 

We begin with the discussion of the rapidity and the transverse momentum distributions 
of the Higgs boson in events with at least one jet, see Fig. 2. The pattern of radiative 
corrections is similar to the hducial cross section case that we just discussed. In the two 
plots in Fig. 2 the relative magnitude of radiative corrections is illustrated in lower panes, 
where ratios of NLO to LO and NNLO to NLO distributions at p = rriH are displayed. We 
will refer to such ratios as iC-factors. We note that similar to the case of the inclusive Higgs 
boson production pp —)■ iL, the NNLO enhancement of the Higgs boson rapidity distribution 
m pp —H + j process is independent of the rapidity. On the contrary, the iF-factors 
for transverse momenta distributions have a more interesting shape. Indeed, we observe 
the instability of da/dp±^H at the value of the Higgs boson transverse momentum equal to 
the value of the jet transverse momentum cut. This is the manifestation of the so called 
Sudakov-shoulder effect [27]. Just above p±^H ~ 30 GeV, the NNLO corrections are small 
but they increase to about 30% at around p±^H ~ 75 GeV and then start to decrease again. 

Next, we consider kinematic distributions of the QGD radiation that accompanies the 
Higgs boson production. The rapidity and the transverse momentum distributions of the 
hardest jet are shown in Fig. 3. Similar to the QGD corrections to the Higgs boson rapidity 
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Figure 4: The distribution of the total transverse energy H± in H + j —?■ 77 production for the 
8 TeV LHC. The inset shows ratios of differential cross sections at different orders in perturbation 
theory for the factorization and the renormalization scales set to the mass of the Higgs boson. 


distribution, the NNLO iC-factor for the hardest jet rapidity distribution is flat. The trans¬ 
verse momentum distribution is re-shaped slightly, with corrections being larger at smaller 
and smaller at larger p±j- The distribution of the total transverse energy H± of QCD 
radiation dehned as shown in Fig. 4. The sum is taken over all jets observed 

i 

in a given events. The NNLO QCD corrections for this observable are smaller for smaller 
iL_L, but, eventually, they increase and flatten out. This is similar to what happens to this 
observable already at next-to-leading order. 

Finally, we turn to distributions that describe kinematic properties of individual photons 
that originate from the Higgs boson decays. In Fig. 5, we show transverse momentum and 
rapidity distributions of the two photons. Similar to other cases, we hnd a uniform enhance¬ 
ment of the rapidity distribution and some shape-dependent NNLO effects in transverse 
momenta distributions. However, the shape-dependence of QCD corrections is signihcantly 
reduced at NNLO compared to the NLO case. This is particularly true for events with trans¬ 
verse momenta at around maximal values for the transverse momenta of both the harder 
and the softer photons. 

In Fig. 6 , we show the distribution of the photon decay angles in the Collins-Soper 
reference frame dehned as 


cosr = I sinh( 7 /.^^ - 1/72)1- 


2p±7]^P±72 


( 4 ) 


+ {p±,H/mHy 

The cos 6* distribution is important for studying the spin-parity quantum numbers of the 
Higgs boson [28]. We hnd that the shape of this distribution is very well predicted by 
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Figure 5: Rapidity and transverse momentum distributions of the harder and softer photons in 
pp ^ H + j at the 8 TeV LHC. The insets show ratios of differential cross sections at different 
orders in perturbation theory for the factorization and the renormalization scales set to the mass 
of the Higgs boson. 

leading order perturbative QCD computations, with both NLO and NNLO QCD corrections 
providing a uniform enhancement. This observation should enable the reduction of the 
uncertainty associated with the modelling of this observable and, perhaps, lead to improved 
limits on exotic features of the observed Higgs resonance. 

Finally, in Fig. 7 we compare the ATLAS measurements with our computations of the 
hducial volume signal in pp —)■ H + j —)■ 77 + j. The inclusive one-jet cross section was 
already discussed at the beginning of this Section; we remind the reader that the result of 
the ATLAS measurement is signihcantly higher than the NNLO QCD prediction for the 
inclusive one-jet cross section. In the left pane of Fig. 7 we present a similar comparison for 
the exclusive jet cross sections. We see that the situation is similar for all jet multiplicities 
and that the discrepancy increases for higher-multiplicity bins. In the right pane of Fig. 7, 
theoretical and experimental results for the transverse momentum distribution of the hardest 
jet are compared. For this observable, the ATLAS results are higher than the theoretical 
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Figure 6 : Distribution of the photon decay angle in the Collins-Soper frame at the 8 TeV LHC. 


The inset shows ratios of differential cross sections at different orders in perturbation theory for 


the factorization and the renormalization scales set to the mass of the Higgs boson. 
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Figure 7: Left pane: comparison of exclusive jet cross sections in pp —)• H j —)• 77 + j computed 
in this paper and measured by the ATLAS collaboration. Right pane: comparison of the leading 
jet transverse momentum distribution. The selection criteria are described in the text. 


predictions in all p_L-bins except one where the experimental error is the largest. It is also 
clear that the shapes of theoretical and experimental distribntions are different. It follows 
from the plots in Fig. 7 that, cnrrently, the ATLAS data is not precise enongh to allow for a 
meaningfnl comparison with available theoretical predictions. This will nndonbtedly change 
once enough luminosity at the 13 TeV LHC is collected. 


B. H ^ VF+fF-^ e+p-nn 


In this subsection, we present the results for the process pp —?■ iL+j —)■ W~^W + j at the 
13 TeV LHC. The selection criteria are described in Section H. In principle, many kinematic 
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Figure 8: Left pane: the production cross section for pp ^ H + j ^ W~^W~ + j —)• e'^pL ~+ j 
at the 13 TeV LHC is shown as a function of the jet transverse momentum cut. The inset shows 
ratios of differential cross sections at different orders in perturbation theory for the factorization 
and the renormalization scales set to the mass of the Higgs boson. Right pane: exclusive jet cross 
sections for pp —?■ + j at the 13 TeV LHC. The selection criteria are described in the 

text. 

distributions that can be studied in the H + j production process are independent of the 
decay mode of the Higgs boson. To avoid overlap with the previous subsection, we present 
here only those distributions that are particular to the W~^W~ hnal state. 

We begin, however, with the discussion of the hducial cross sections. We hnd 

(Tlo = fb, n^Lo = fb, ctnnlo ~ 21.9l5'7 fb. (5) 

In general, the perturbative expansion of the 13 TeV cross sections is similar to what was 
observed at 8 TeV. At p = ttih, the NLO cross section is larger than the LO cross section by 
43% and the NNLO cross section exceeds the NLO cross section by 18%. In Fig. 8 we display 
results for cross sections as a function of the jet transverse momentum cut and exclusive 
jet cross sections at different orders in perturbation theory. The behavior of the exclusive 
one-jet cross section at 13 TeV is slightly worse than at 8 TeV. We attribute this to a too 
small scale variation at NLO, which leads to the NNLO result for the one-jet cross section 
being outside of the NLO scale variation band. 

Selection criteria for the Higgs signal in if —)■ W^W~ as well as analysis of anomalous cou¬ 
plings in this process rely on kinematic distributions of charged leptons. A good understand¬ 
ing of these distributions is therefore important. In Fig. 9 we show the transverse momentum 
distribution of a positively charged lepton in pp —)■ iL +j —>■ W^W~ +j —>■ e~^p~u9+j at the 
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Figure 9: Left pane: the transverse momentum distribution of a positively charged lepton in 
pp —>■ + j at the 13 TeV LHC. Right pane: the distribution of the azimuthal opening 

angle of the two leptons in pp —)■ H{e'^p~vi') + j at the 13 TeV LHC. The selection criteria are 
described in the text. The insets show ratios of differential cross sections at different orders in 
perturbation theory for the factorization and the renormalization scales set to the mass of the 
Higgs boson. 

13 TeV LHC and the azimuthal opening angle distribution of the two leptons. In both cases, 
QCD radiative corrections do not change the shapes of the distributions significantly. The 
distribution of the invariant masses of the two leptons rrti+i- and the transverse mass m±_ 
are displayed in Fig. 10; the NNLO QCD corrections to those distributions are remarkably 
uniform. 

IV. CONCLUSIONS 

We extended the recent NNLO QCD computation of the H + j production in proton 
collisions by including decays of the Higgs bosons to electroweak gauge bosons iL —)■ 77 , 
H —)■ W^W~ and H — ZZ. Leptonic decays of Z’s and LF’s, with all spin correlations, 
are fully accounted for. This allows us to calculate fiducial volume cross sections and vari¬ 
ous kinematic distributions through NNLO in perturbative QCD in a manner that is fully 
consistent with selection criteria applied in experiments. In particular, it becomes possible 
- for the hrst time - to confront hducial volume studies of the pp ^ H + j —>■ 77 + j process 
performed by the ATLAS collaboration at the 8 TeV LHC [5] with NNLO QCD predictions. 

We presented a number of results for fiducial volume cross sections, acceptances and vari¬ 
ous kinematic distributions for both inclusive and exclusive if ( 77 ) +j production processes. 
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Figure 10: Left pane: the lepton invariant mass distribution in pp —)• H{e'^p~i'u) + j at the 
13 TeV LHC. Right pane: the W~^W~ boson transverse mass distribution in pp —)> H{e~^ui') +j 
at the 13 TeV LHC. The selection criteria are described in the text. The insets show ratios of 
differential cross sections at different orders in perturbation theory for the factorization and the 
renormalization scales set to the mass of the Higgs boson. 


We find no indication that pertnrbative QCD breaks down and reqnires resnmmation for the 
jet cut as low as 30 GeV. We have also studied the WW hnal state in H + j production at 
the 13 TeV LHC. We hnd that most of the kinematic distributions used to distinguish this 
channel from backgrounds, show uniform enhancement when NNLO QCD corrections are 
included. Changes of shapes of such distributions - if any - are already properly captured 
by the NLO QCD computations. 

As a hnal remark, we note that the availability of higher order QCD predictions for 
hducial volume quantities should allow direct and precise studies of the ratios of Higgs 
signals. The idea that ratios of cross sections are useful for reducing theoretical uncertainties 
is, of course, well-known and appreciated. However, given the availability of the NNLO 
QCD computations for hducial cross sections, no extrapolations should be required. As an 
illustration we compute the ratio of hducial cross sections for pp —)■ iL j —>■ 77 -|- j at the 8 
TeV LHC and the pp —)■ iL -|- j —)■ WW* + j ^ e+/i“izp -f j at the 13 TeV LHC. We obtain 


R 


a 


wwh-i 


WW^e+fi uu,l 3 TeV 
H+j 


_ O OQ —0.06 Q 00 —U.U4 o 00 - 

— Z.Oa^gQ 4 , Z.OO_,_Qg 5 , Z.OZ_,_gQ 2 , 


>-0.04 


1-0.04 


a 


77,8 TeV 
H+j 


( 6 ) 


at leading, next-to-leading and next-to-next-to-leading order in perturbative QCD, respec¬ 
tively. The convergence of the series is striking; at NNLO QCD, we are able to predict 
Rwwi'j'y with the precision of better than 2 %. Since Rwwf-y'y is proportional to the ratio 
of the Higgs couplings to two photons and to W-bosons, confronting a precise prediction 
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for this observable with results of experimental measurements should allow for stringent 
constraints on the deviations of these couplings from their Standard Model values. 
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